On properties of elementary excitations in fractal media 
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It is shown that elementary excitations in fractal media obey the so-called parastatistics of a variable 
order. We show that the order of the parastatistics N{k) is a function of wave numbers k which 
depends on the fractal dimension D as N{k) ~ k^~^ and represents a specific characteristic of such 
media. This function N{k) defines properties of the ground state for excitations and the behavior 
of the spectrum of thermal fiuctuations. In particular, in fractal media fiuctuations of the density 



acquire an amplification by the factor N{lj) 



which can be related to the origin of l//-noise. 



It is well known that a physical system in which parti- 
cle distribution follows, in some range of scales, a fractal 
law N (r) ^ (where N (r) is the number of parti- 
cles contained within a radius r) with 13 < 3 is not an 
exotic example, rather it represents a typical situation 
(e.g., see Refs. [1] - [3] and references therein). Neverthe- 
less, an acceptable quantitative theory which allows for 
the description of thermodynamical properties of such 
systems is still absent. In the present paper we, however, 
show that the possibility to construct such a theory has 
recently appeared. First, we note that thermodynami- 
cal properties of a system are determined by properties 
and behaviour of elementary excitations in the system. 
In quantum theory excitations are described in terms of 
respective particles, while the number of degrees of free- 
dom of a particular excitation can be characterized, in 
the Fourier space, by the spectral number of modes in 
the interval of wave numbers between k and k + dk as 



(1) 



The behaviour of the function v (k) depends on the di- 
mension of the system as ly' (k) ^ k^ and for regular (or 
homogeneous) systems the exponent D = n takes inte- 
gral values n = 1,2, or 3. In the case of fractals, how- 
ever, the dimension D is not an integral number, while 
in a more general case (e.g., in the case of multifractals) 
this function may change its behaviour for different in- 
tervals of scales, i.e., D = D{k). Moreover, a change 
in the thermodynamical state of a system may lead to 
essential transformations in the behaviour of the func- 
tion {k) (e.g., in the case of percolation systems [2]). 
Thus, an adequate description of such systems requires 
that the spectral number of modes ly (fc) has to be a vari- 
able which, in general, depends on wave numbers and, in 
the presence of self-organization processes, on time. 

In the standard field theory the spectral number of 
modes is always fixed by the dimension of space. How- 
ever, there exists a modification of the field theory 



(MOFT) in which the number of fields (and, therefore, 
of field modes) can vary. Such a modification was pro- 
posed in Ref. [4] to account for spacetime foam effects 
and it was recently developed in Refs. [5]. In particular, 
MOFT was shown to predict a significant amplification 
of the gravitational force at galaxy scales (which is really 
observed and represents the so-called dark matter prob- 
lem) and a fractal law for the galaxy distribution with 
dimension Z3 « 2 [5] which is in a full agreement with 
the observed picture of our Universe [6]. Thus, we can 
expect that the theoretical scheme of MOFT provides a 
natural basis to describe excitations in systems with the 
fractal behaviour. 

The measure (1) can be used to define the number 
of field modes N (fc) as v (fc) = k^N (k) / {4tt^) (to be 
specific we consider a fractal volume). In the standard 
field theory TV (fc) = 1, while in a general case N (k) 
is a variable which depends on the wave number k and 
time (in particular, the fractal law gives N (fc) ~ k^^^). 
From the mathematical standpoint this means that exci- 
tations obey the so-called parastatistics which generalizes 
the standard Bose and/or Fermi statistics and the num- 
ber of field modes N {k) plays the role of the order of the 
parastatistics. We note that the first attempt to gener- 
alize the Bose and Fermi statistics was made by Gentile 
[7], while the parastatistics was suggested by Green in 
Ref. [8]. Since then, the parastatistics has been studied 
in many papers (e.g., see Refs. [9] - [11]). In particu- 
lar, experimental observation of the fractional quantum 
Hall effect [12] stimulated the interest in various nontra- 
ditional statistics [13,14] and an attempt was made to 
relate the parastatistics to high temperature supercon- 
ductivity [15]. However, the parastatistics is still con- 
sidered as an exotic possibility and has not received an 
adequate attention (e.g., see Ref. [16]). MOFT [5] gives 
a specific realization of the parafield theory in which the 
order of the parastatistics represents an additional vari- 
able. The discussion above shows that the parastatistics 
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of a variable order has a straightforward relation to the 
formation of properties of elementary excitations in frac- 
tal media. 

Consider as an example sound waves in a system of 
identical g— particles whose distribution in space has, at 
low temperatures T — > 0, an irregular (i.e., fractal) char- 
acter po (x) . We note that the case of fcrmionic ex- 
citations can be considered in the same way. Sound 
waves represent oscillations in the density 6p{x,t) = 
p {x, t)~'Po (x) and the velocity Su {x, i) = u {x, t) — UQ (x) 
of particles (where po, uq are mean values which, in frac- 
tal media, have a specific dependence on coordinates). 
In what follows we denote these functions as a field 
A {x, t) = {dp {x, t) , 6u {x, t)). This field can be expanded 
in modes 

A{x,t) = ^Yl "-^."e" (^) e^"' + <a< (k) e-'=^ (2) 

k,a 

where V is the volume and Ca (fc) is the polarization vec- 
tor (we assume that any additional normalization factor 
is included in Ca {kj). The coefficients a^^ and a^^a play 
the role of creation and annihilation operators for ele- 
mentary excitations (phonons). Thus, in the standard 
picture these operators obey the relations 

ak,aak',l3 ~ ^t',0^k,oi = dk,k'da,0, (3) 

which mean that phonons represent Bose particles, while 
the Hamiltonian for free particles is given by 

Ho = ^Wfc,c«afc,aafe,a , (4) 

k,a 

where u)k,a is the energy of a phonon. 

As it was discussed above in fractal media oscillations 
have a more complex character than in homogeneous me- 
dia (e.g., like in ideal crystals or gases) and for an ade- 
quate description of oscillations a single field A (x, t) is 
not enough. Indeed, to illustrate this statement we can 
imagine the situation when the system of particles splits 
into a set of almost independent thin surfaces (e.g., the 
system represents a fractal with dimension D ~ 2), so 
that excitations for every surface can be described by 
its own field Ai {x,t) (where i numerates the surfaces). 
Such surfaces can have a chaotic distribution in space 
and, therefore, to distinguish them in an explicit way is 
not so trivial. In real fractal media such a splitting works 
in the Fourier space (e.g., see Ref. [5]) and, therefore, in 
a general situation we can only state that particles of 
a system can be divided in groups n = {ni,n2, ■■■Un} 
{J2i — where n is the total munber of particles), 
so that every group is responsible for the formation of its 
own excitations. In general, such a decomposition can 
depend on time, while its character and the distribution 
of the groups in the space of modes A'' (fc) constitute par- 
ticular properties of the medium which require the exact 



knowledge of the dynamics of particles (in simple cases, 
however, this information can be found from simple ther- 
modynamical considerations, e.g., see Ref. [5]). As we 
will see, the distribution of groups N (fc) in the space 
of modes defines spectrum of thermal fluctuations in the 
system and, therefore, for practical needs this function 
can be determined from direct measurements . 

It is important that g— particles are identical, and, 
therefore, the decomposition n = ^ rij has a conditional 
character. In particular, the mean values po and uq are 
determined by all groups, while perturbations acquire 
analogous decomposition A = X^Aj. The identity of 
particles q results in the analogous identity of fields Ai. 
Thus, from the mathematical standpoint excitations will 
be particles which obey the parastatistics, while the num- 
ber of fields (groups) N plays the role of the order of the 
parastatistics. 

In what follows we, for the sake of simplicity, neglect 
the presence of the index a (which numerates different 
polarizations). In this manner, when the number of fields 
is variable, the set of creation /annihilation operators 
{ak, a^} is replaced by the expanded set {at (j) , (j)}, 
where j S [1, ...N (A:)], and N (k) is the total number of 
fields (groups) for a given wave number k. Then the 
energy for a free field can be written as 

N(k) 

^O^J^Jl '^^"'k 0) O'k (j) • (5) 
k j=l 

The fields arc; identical and are supposed to obey the 
same statistics as that of g— particles. It is conve- 
nient to consider the set of creation/annihilation oper- 
ators for field modes {C+ (n, k) ,C {n,k)}, where k is the 
wave number and n is the number of particles in the 
given mode. Since the number of modes corresponds to 
the true degrees of freedom there exists a relation be- 
tween the operators C"*" (n, k) and C (n, k) and the stan- 
dard creation/annihilation operators tp'^ (x) and ip (x) for 
g— particles, i.e., C{n,k) = f A{n,k,x)ip {x) (Px with 
some kernel A{n,k,x). Thus, operators C+(n, fc) and 
C (n, k) obey the relations 

C (n, k) C+ (m, k') ± C+ (m, k') C (n, k) = Sn„,Skk' , (6) 

where the sign reflects the statistics of g— particles. We 
note that in the case of the so-called parafield theory 
there exists a spin-statistics theorem [If] (analogous to 
the Pauli theorem) which states that bosonic modes 
should be quantized according to the Fermi statistics, 
while fermionic modes should obey the Bose statistics. 
This means that the fractal behaviour can be associ- 
ated only with fermions. In the case of a fractal medium 
in which atoms have an integral spin we should assume 
that the distribution of modes is determined by electrons. 
Mathematically, this can be expressed by the fact that in 
the case when g— particles are bosons the necessary ker- 
nel A (n, k, x) does not exist. Thus, in what follows we 
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consider the case of the Fermi statistics for modes only 
(the sign + in (6)). 

The eigenvalues of the Hamiltonian (5) can be written 
straightforwardly 



Ho = ojknN (n, k) . 



(7) 



where N (n, k) = C"*" (n, k) C (n, k) is the number of field 
modes with the given wave number k and the number of 
phonons n. The fact that modes obey the identity prin- 
ciple results in a redefinition of the standard creation 
and annihilation operators for phonons (see, for details, 
Ref. [4]). Indeed, since fields Ai (corresponding to dif- 
ferent groups of g— particles) are identical they cannot 
be detected separately. What is really measured in ex- 
periments is the total perturbation A = Y^Ai which is 
given by the expression (2) in which, however, we should 
replace the creation/ annihilation operators ak and by 



sums 



N{k) 



N{k) 

E 4 (i) • 



(8) 



In the representation of occupation numbers N (n, k) 
these operators have the expression [4] 

ak = YVrr+lC+ {n,k)C{n+l,k), (9) 



a+ = J2 Vn + TC+ (n +l,k)C {n, k) 



Thus, from (6) we find that these operators obey the 
relations 



ttka^, - a^ittk = N {k) 5k,k' , 



(10) 



which generalize the standard relations (3) and are a con- 
sequence of the fact that phonons obey the parastatistics. 
The order of the parastatistics A'' (fc) is the operator of 
the total number of field modes for a given wave number 
k which is given by 



N{k) = J2^iri,k). 



(11) 



Physically, this operator characterizes the distribution 
of the number of groups of g— particles in the space of 
modes, while A'' = ^ i^) gives the total number of 
particles in the system. 

In the thermal equilibrium state (and in the absence of 
self-organization processes) the number of fields is con- 
served and the operator N (k) can be considered as an 
ordinary function of wave numbers. Thus, we assume 
that N {k) represents a structural function of the system 
which should remain constant (at least, for not very high 



temperatures) . Then the thermodynamically equilibrium 
state will be characterized by mean values for occupation 
numbers of the type 



(7V(fc,n))=(exp(l^^^^)+l 



(12) 



where the chemical potential fik can be expressed via the 
structural function N {k) as 



N{k)=Y,U^p( 

n. ^ ^ 



nujk — Hk 



+ 1 



(13) 



In the limit T — > the system reaches the ground state. 
From (12) we find that the ground state is characterized 
by the occupation numbers 



N (n, k)=9 {nk - nojk) 



(14) 



where 9 (x) is the Heaviside step function. Thus, in the 
limit T — » we find that the chemical potential /x/j is 
expressed via the structural function N (k) as 



N{k) = {nk - nwk) = 1 + 



(15) 



where [x] denotes the integral part of the number x. 

From (14) we find that the ground state is character- 
ized by a nonvanishing number of phonons no (fc) 



oo ^ 

no (fc) = E fc) = 9^ (fc) (^) 



1) 



(16) 



n=0 



and corresponds to the energy = J2^kno (k). This 
accounts for the fact that in the limit T ^ the distribu- 
tion of g— particles in the medium has an inhomogeneous 
(fractal) character (i.e., in the ground state the density of 
(7— particles contain a specific dependence on coordinates 
po {x), e.g. see Ref. [5]). 

The expression (15) allows to understand some gen- 
eral properties of the structural function N{k). Indeed, 
consider the case when fik = M- Then there exists a re- 
gion of wave numbers > /i in which N (k) = 1 and 
phonons behave as standard Bose particles. Thus, the 
wave number k* (at which lu = 11) characterizes the scale 
on which effects of the parastatistics start to show up. In 
this sense regular systems can be considered as those sys- 
tems in which this scale £* = 2'n/k* exceeds the maximal 
possible wavelength. At scales k <^ k* (i.e., as tOk 0) 
we get N (k) ~ /u/wfe which gives the fractal dimension 
£) « 2. 

Now we examine the problem how the structural func- 
tion N {k) can be expressed via the correlation functions 
and hence how it can be directly measured. Substituting 
the expressions (9) for the creation/annihilation opera- 
tors in (2) we find that correlation functions have the 
structure 
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(^W,A(. + r)) = ^^°°$^W^f (17) 
where 

^^k) = k^\e{k)\''{2nk + N{k)) (18) 
and n/s = (a^ak) is given by 

nfc = ^ (n + 1) {N (n + 1, k)) (1 - {N (n, A;))) , (19) 

n 

where the occupation numbers (TV (n, k)) arc defined in 
(12). For small temperatures T ^ we get ^ 
and the correlation function (fc) explicitly defines the 
structural fimction of the system N (k). 

The explicit form of the function rik can be found 
as follows. We recall that we assume that the num- 
ber of modes N (fc) represents a constant function (the 
structural function of the system). Then from (10) 
we see that the difference between operators and 
ak and the standard creation and annihilation opera- 
tors for Bose particles comes from the coefficient Uk = 
y/N {k)ck (where Cfe is the standard annihilation oper- 
ator, i.e., [cfe.c^,] = Sk.k')- Hence for thermal equi- 
librium state we should get the expression {a^ak) = 
N (k) (c+Ck) = N (k) (exp (^) - iy\ Therefore, from 
the formal standpoint the fact that phonons obey the 
parastatistics can be accounted for by the substitution 
{2nk + 1) = coth (^) N (k) coth (^) in all thermo- 
dynamical formulas and, in general, the structural func- 
tion N (k) defines deviations of thermodynamical poten- 
tials and fluctuations from the standard thermodynam- 
ical laws. In particular, for phonons we get uj — ku, 
where u is the sound speed (in general u can depend on 
the polarization of phonons, but we neglect here this ef- 
fect). Thus, if in some range of scales a system possesses 
a fractal behaviour with dimension D, wc find that in 
that range fluctuations will acquire an amplification by 
the factor N (oj) ~ w"" with a = 3 — D. Such kind of 
amplification is indeed observed in various systems and 
represents the so-called 1/f -noise problem [17]. Hence 
we may conclude that at list in some class of systems the 
origin of 1// —noise can have straightforward relation to 
the fractal behaviour and, therefore, to the parastatis- 
tics of elementary excitations. This problem, however, 
requires the further investigation. 

In this manner, wc have shown that in a wide class of 
systems (which demonstrate an anomalous behavior for 
thermodynamical potentials and fluctuations, presum- 
ably at low frequencies) elementary excitations obey the 
so-called parastatistics of a variable order N{k). The 
order of the parastatistics N (fc) reflects properties of a 
system at low temperatures (e.g., the inhomogeneous dis- 
tribution of particles) and represents a specific, for every 
given system, structural function which can be deter- 
mined by means of direct measurements of fluctuations. 



In fractal media this function behaves as N (k) ^ k^^^, 
while for a more general case (i.e., in general disordered 
systems) this function can have an arbitrary behaviour. 
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